Testing Multi-Field Inflation: A Geometric Approach 



Courtney M. Peterson^, Max Tegmark^ 
^ Dept. of Physics, Harvard University, Cambridge, MA 02138, USA 
^ Dept. of Physics & MIT Kavli Institute, Massachusetts Institute of Technology, Cambridge, MA 02139 

(Dated: November 3, 2011) 

We develop an approach for testing inflation models with multiple scalar fields by linking geomet- 
ric and kinematical features of their inflationary Lagrangians to observable quantities like the power 
spectra, bispectrum, and trispectrum. Our approach also provides geometric intuition for when a 
complicated multi-field model can be well- approximated by a model with one, two, or a handful of 
fields. To arrive at these results, we focus on the mode interactions, simplify them using a novel re- 
sult, and then explore how these interactions depend on the geometry of the inflationary Lagrangian 
and on the kinematics of the associated field trajectory In the process, we introduce a multi-field 
observable {32 that can potentially distinguish two-field scenarios from scenarios involving three or 
more effective fields. We also present a multi- field consistency relation, which involves the primor- 
dial bispectrum parameter /atl, trispectrum parameter tatl, and other spectral observables. These 
combined results provide better intuition into how features in multi-field inflationary Lagrangians 
translate into cosmic observables. 



I. INTRODUCTION 

Inflation solves cosmic conundrums such as the hori- 
zon, flatness, and relic problems [iHS]. It also offers a 
mechanism for producing the primordial density fluctua- 
tions. According to the inflationary paradigm, our Uni- 
verse experienced an early period of quasi-exponential 
expansion. As a result of this inflationary expansion, 
quantum fluctuations were stretched beyond the causal 
horizon. Once beyond the horizon, they became locked 
in as classical perturbations, which eventually initiated 
the formation of galaxies and large-scale structure [6]- 
[TT] . Generically, inflation predicts that these classi- 
cal perturbations should be reflected in a nearly scale- 
invariant spectrum of primordial density fluctuations at 
early times. 

So far, measurements of the Cosmic Microwave Back- 
ground (CMB), large-scale structure, supernovae, and 
gravitational lensing support the inflationary paradigm. 
These measurements reveal a nearly scale-invariant spec- 
trum of small primordial fluctuations, the presence of 
superhorizon fluctuations, and that our Universe is es- 
sentially flat, as predicted by inflation (see [12] and ref- 
erences therein). But ultimately, we want to use cosmic 
data not just to test the inflationary paradigm, but to 
find the particular inflationary model that describes our 
Universe. 

Of the myriad inflationary models, there is good rea- 
son to consider models where inflation is driven by mul- 
tiple scalar fields. Many theories beyond the Standard 
Model — such as grand unification, supersymmetry, and 
effective supergravity from string theory — predict the ex- 
istence of multiple scalar fields, which makes the presence 
of multiple fields likely during the hot, early Universe. 
Moreover, multi- field models have become increasingly 
popular in recent years. These two facts strongly moti- 
vate the study of multi-field inflation. 

But the sobering reality is that there is a staggeringly 
large number of multi-field models, making it impracti- 




FIG. 1. Examples of three-field inflation trajectories that can 
be accurately approximated by an effective model with (a) 
one, (b) two, and (c) three fields, respectively For example, 
the field trajectory in (b) requires more than one effective 
field to represent the trajectory because it curves, but the 
trajectory curves only in a plane, so only two effective fields 
are needed. 



cal to test each model against cosmic data. Indeed, we 
illustrated this point in [13 by examining a few classes 
of models (for example, double quadratic models) and 
rigorously testing them by considering the full range of 
initial conditions. For each type of model, we tested 
more than 10,000 scenarios by varying both a parame- 
ter value in the Lagrangian and the initial conditions, in 
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order to constrain that class of models using WMAP con- 
straints on the power spectra. But as we demonstrated 
in [13], such an approach is extremely time-consuming. 
Instead of testing each inflationary scenario one by one, a 
more fruitful approach is to determine how key features of 
multi-field infiationary Lagrangians translate into observ- 
able quantities. For example, how is the geometry of the 
infiationary potential refiected in the cosmic observables? 
Is there a way to tell from cosmic data whether a one- 
field or two-field model can fit all measurements, as illus- 
trated in Figure 1, or whether more fields are required? 
And what is the role of nonstandard kinetic terms in 
determining the cosmic observables? In this paper, we 
attempt to answer these and other questions, with our 
ultimate goal being to understand how cosmic data can 
constrain the form of the infiationary Lagrangian. 

We do this by building on the pioneering work of oth- 
ers. Initial work on understanding the perturbations and 
power spectra generated by general multi-field infiation 
was done in p!4H28] . The specific case of two-field infia- 
tion was treated in [TSj [28ti4Q] . Work towards calculating 
other spectra, such as the bispectrum and trispectrum, 
in general two-field and multi-field infiation was done by 
[4r-^62^, among others. Interestingly, while developing 
generic semianalytic formulas for the spectra from multi- 
field infiation has received much attention, the sourcing 
relations among the modes in multi-field models has re- 
ceived less attention. Understanding the sourcing effects 
is important as they provide the link between the infia- 
tionary Lagrangian and the cosmic observables. This lack 
of attention is at odds with investigations of two-field sce- 
narios, in which the sourcing effects have been treated in 
detail. Related to the mode sourcing is the effective num- 
ber of fields needed to describe an infiationary scenario. 
Knowing the effective number of fields is important be- 
cause it can lead to dramatic simplifications of the ex- 
pressions for the spectra and it allows us to focus on the 
most important features of an infiationary model. For 
example, if a multi-field model with ten scalar fields be- 
haves like a two-field model during the last several dozen 
e-folds of infiation, then using this fact greatly enhances 
our understanding of that model, as well as largely re- 
duces the computational complexity of the spectra. 

In this paper, we focus on the underappreciated series 
of mode sourcing relations and how they are determined 
by the infiationary Lagrangian, in order to extend our 
insights from [13, 69 to multi-field infiation. The rest of 
this paper is organized as follows. In Section [llj we cover 
the dynamics and kinematics of the background fields, 
and we discuss slow-roll and other limiting cases of be- 



havior. In Section III, we present equations of motion for 
the field perturbations in both the given and kinematical 
bases. We then prove that most of the coefficients of the 
mass matrix vanish in the kinematical basis, and we use 
this fact to derive a simplified series of sourcing equations 
for the modes. In tandem, we discuss how the geometric 
and kinematical features of the infiationary Lagrangian 
determine the interactions among the modes. Finally, 



in Section [TVj we use these sourcing equations to exam- 
ine the effective number of fields in multi-field models 
and to explore how this number is refiected in spectral 
observables. We also generalize our two-field semiana- 
lytic formulas for the bispectrum and trispectrum [69^ to 
multi-field infiation, introduce a new multi-field consis- 
tency condition, and identify a spectral observable that 
can be used to distinguish two-field models from models 
with three or more fields. This work paves the way to- 
ward a better understanding of how the features of multi- 
field infiationary Lagrangians are refiected in the cosmic 
observables. 



II. BACKGROUND FIELDS 

We start our treatment of multi-field infiation by ex- 
amining the dynamics and kinematics of the background 
fields. In Section [TTAI we review notation and the equa- 
tion of motion for the background fields. In Section II B[ 
we outline a framework for understanding the field vector 
kinematics, which is based on work by [13^ 21, 23 , 30 . We 
conclude our treatment of the background fields in Sec- 
tion II C| by exploring the standard slow-roll limit and 
other limiting cases of behavior. This extended work- 
up — particularly on the background kinematics — lays the 
necessary groundwork for us to connect the geometric 
features of an infiationary Lagrangian to the cosmic ob- 
servables. 



A. Evolution Equation for the Fields 

We consider infiationary scenarios driven by an arbi- 
trary number of scalar fields, 0^, where i = 1,2, 
and d is the total number of scalar fields present during 
infiation. Here, we allow for the possibility that some 
number of the d scalar fields do not contribute to the in- 
fiationary expansion. Which of these fields do and do not 
contribute to the infiationary expansion does not need to 
be known in advance, as the formalism presented here 
will automatically take care of this. To represent the 
fields compactly, we use the notation 



= (^1,^2, 



(1) 



and we call the field vector for short, even though the 
fields do not transform as vectors. Here, we express the 
fields in units of the reduced Planck mass, ffi = 
and we work in units where c = h = fh = 1. 

Assuming Einstein gravity, we take the non- 
gravitational part of the infiationary action to be 



2^ ^''^x^^^x- 



-gd^x, (2) 



where (/^^ is the spacetime metric and Gij determines the 
form of the kinetic terms in the Lagrangian. G,j is called 
the field metric, which we assume is a function of only the 
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fields, and it can be viewed as inducing a field manifold. 
If the kinetic terms are canonical, then Gij = Sij^ and 
the field manifold reduces to Euclidean space. Here, we 
allow both the field metric and the inflationary potential, 
V{(f)), to be arbitrary. 

Before proceeding, we introduce some notational short- 
hand. We use Latin indices to represent quantities re- 
lated to the fields, For vectorial quantities lying in 
the tangent and cotangent bundles of the field manifold, 
we use boldface vector notation and standard inner prod- 
uct notation. The inner product of two vectors A and B 



where F*^ and Vj are the Levi-Civita connection and 
the covariant derivative, respectively, associated with the 
field metric. Therefore, 77 represents the covariant rate of 
change of the field velocity vector with respect to the field 
manifold — that is, it represents deviations from perfect 
parallel transport of cj)' . By working in terms of D and 
the covariant derivative V, we are able to write all the 
equations of motion in manifestly covariant form. 



Field Vector Kinematics 



A B 



(3) 



where we use the symbol ' on a naturally contravariant 
or covariant vector to denote its dual, e.g., cj)^ = {Gijcj)^) 
and = {G'^Vj). Norms are defined as the square root 
of the inner product of a vector with itself: 



VAtA. 



(4) 



Lastly, instead of working in terms of the coordinate time, 
t, we work in terms of N ^ which represents the logarith- 
mic growth of the scale factor, a{t)\ 



dN = d\na = H dt, 



(5) 



where H = ^ is the Hubble parameter. In other words, 
represents the number of e-folds of the scale factor. We 
work in terms of N because it is dimensionless, it relates 
to a more physical measure of time, and it simplifies the 
equations of motion [131 IE]- To denote differentiation 
with respect to we use the notation 



_d_ 
dN' 



(6) 



We are now ready to present an equation of motion 
for the fields — the main goal of this section. By impos- 
ing covariant conservation of energy, we derived such an 
equation using as the time variable in [13]. We found 
that the equation of motion for the fields can be written 
as 



V 



(3-6) 



where 



e = -(lnif)' = i0'.0', 

and the covariant field acceleration 77 is defined as 



V 



dN 



(7) 



(8) 



(9) 



The symbol D acting on a contravariant vector X' means 



(10) 



The kinematics of the background fields are important 
because they allow us to follow how the geometric fea- 
tures of an inflationary Lagrangian translate into values 
for the spectral observables. Conversely, experimental 
bounds on spectral observables often translate into con- 
straints on the background kinematics, which in turn con- 
strain the inflationary Lagrangian. In other words, the 
field kinematics serve as intermediaries between the infla- 
tionary Lagrangian and the spectral observables. We dis- 
cussed this in great detail for the case of general two-field 
inflation with noncanonical kinetic terms in [13 . Here, 
we extend that and other work by presenting a kinemati- 
cal framework for the general case of multi-field inflation. 

The first kinematical framework for multi-field infla- 
tion was constructed by [30] for the specific case of two- 
field inflationary potentials with canonical kinetic terms. 
Soon after, Nibbelink and Van Tent extended this frame- 
work to describe multi-field inflation with noncanonical 
kinetic terms [21j i23^ . In [13^ , we slightly modified this 
framework in applying it to two-field inflation with non- 
canonical kinetic terms. Here, we combine elements from 
[131 [211 113 in presenting a kinematical framework for 
multi-field inflation. We will note where any quantities 
we introduce are different from those originally proposed 
by Nibbelink and Van Tent. 

To construct a kinematical framework, we start in the 
normal manner by considering the velocity, acceleration, 
and higher-order derivatives of the coordinates. Here, 
the coordinates are the scalar fields, which represent the 
coordinate position on the manifold induced by the field 
metric. In analogy to Newtonian mechanics, cj) represents 
the position, 0' is the velocity, and 77 = represents 
the covariant acceleration, where ^ is effectively defined 
by Eq. (10). Similarly, we can define higher-order covari- 



ant derivatives of the field velocity ad infinitum. The jerk 
is defined as 



dm ' 



(11) 



And an equation of motion for the jerk can be obtained 
by differentiating Eq. ([7| once, which yields 



(3-6) 



M 



(3-e)2J 



(12) 
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where the mass matrix, is defined as 



(13) 



and is symmetric. Similarly, we represent the (n — l)-th 
covariant derivative of the velocity by the notation 



An) 



D 



dN J 



(n-l) 



(14) 



and an equation of motion for x'"^'' can be obtained by 
differentiating Eq. ^ a total of n — 2 times. For com- 
parison, Nibbelink and Van Tent [2^^ defined a series of 
higher-order kinematical vectors as 



2)("-i)0; 



(15) 



where ; represents the derivative with respect to the arbi- 
trary time variable r and V is the "slow-roll derivative" Q 
Both constructs have their utility: Nibbelink and Van 
Tent's construct makes their vectors manifestly indepen- 
dent of the choice of time coordinate, while our construct 
is more physically intuitive and can be used to simplify 
many infiationary expressions to a greater degree. 

To accompany these kinematical vectors, one can con- 
struct a basis that is induced by the field vector kinemat- 
ics. Working in this basis simplifies the equations of mo- 
tion for the field perturbations. In [13^, we dubbed this 
kinematically-induced basis the kinematical basis. Gor- 
don et al. \30' first developed such a basis for two-field 
infiation with canonical kinetic terms, and Nibbelink and 
Van Tent extended this basis to multi-field infiation with 
noncanonical kinetic terms ^21||23j. In [13^, we adopted 
this same basis, but used N as the time variable when we 
applied it to two-field infiation with noncanonical kinetic 
terms. The construction of this basis is as follows. The 
first basis vector, ei, is chosen to lie in the direction of the 
field velocity, so that it is parallel to the field trajectory. 
The second basis vector, e2, is constructed to lie along 
the part of the field acceleration that is orthogonal to the 
field velocity, in the direction that makes e2 • ?7 > 0. Us- 
ing the Gram-Schimdt orthogonalization procedure, this 
process is continued until there are d basis vectors. In 



The "slow-roll derivative" is defined as V(b'^A) 



y dT 



where b = —goo and A is independent of b. Our 
kinematical vectors differ from Nibbelink and Van Tent's in two 
ways: (1) the effective order in the slow-roll expansion, which 
differs because of the factor of |</)'| in the denominator in Eq. 
( |15| , and (2) the expressions themselves — that is, our series dif- 
fers from theirs even when the order of Eq. ( |15| is adjusted by 
multiplying by |0' |. For more information, we refer the interested 
reader to 1231. 




FIG. 2. An example showing how the kinematical basis is 
constructed from the kinematics of the background fields. The 
green curved path represents the trajectory of the field vector 
in a three-field inflation scenario, as a function of the time 
variable N. At each point on the trajectory, the ei basis 
vector is chosen to point in the direction of the field velocity, 
62 points in the direction of the perpendicular acceleration, 
and es is constructed to be orthogonal to the first two basis 
vectors. 



other words, the kinematical basis vectors are 



ei 



e2 



(I-eie|)r/ 
^I(l-eiel)r7| 



(16) 



-1 



-1 

i=l 



where I is the appropriate identity matrix. Now, if the 
kinematical vector x*^^^ already lies in the subspace de- 
fined by the basis vectors ei, e2, e^-i, then it is not 
possible to find a projection of that represents a 
new direction in field space. In this case, can simply 
be constructed at will so that it represents a new direc- 
tion that is orthogonal to the subspace spanned by the 
basis vectors ei through e^-i, and then the orthogonal- 
ization process can naturally proceed again. Our basis 
vectors are effectively identical to those of Nibbelink and 
Van Tent [21^ 23 , despite our different definition for the 
kinematical vectors. 

With these basis vectors, we can take projections of 
vectors and matrices. For example. 



An — X 



(17) 
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represents the projection of the kinematical vector 
onto the basis vector e^. Note that because of the def- 



inition of the kinematical basis vectors in Eq. (16), if 



n > m, then Xn^^ = 0- That is, in the kinematical basis, 
(p^ has the sole nonzero component 

v=\cf>'\; (18) 

T] has nonzero components r]i and 7^2; C has nonzero com- 
ponents ^1, ^2, and ^3; and so on. The projection of any 
vector along ei is particularly noteworthy, as it represents 
the vector component parallel to the field trajectory. To 
extend these projections to any matrix A, we use the 
shorthand notation 



for the matrix coefficients and 



= (l-eiel) A(l-eiel) 



(19) 



(20) 



for the above special matrix projection. 

The last useful part of this kinematical framework is 
the derivatives of the basis vectors with respect to the 
number of e-folds of inflation, N. These temporal deriva- 
tives represent how quickly the basis vectors are covari- 
antly changing direction with respect to the field mani- 
fold. In particular, the derivative of the basis vector par- 
allel to the field velocity, ei, represents how quickly the 
field trajectory itself is covariantly changing direction. It 
is given by differentiating Eq. (16), which yields 



Dei m 



-e2. 



dN V 

Similarly, differentiating the second basis vector gives 



(21) 



De2 6 V2 
dJy r]2 V 

The derivative of the nth basis vector is 



Den _ Xn+1 



dN 



in) 
Xn 



X^ 



(n) 



(n-1) 
Xn-1 



(22) 



(23) 



Thus, another benefit of the kinematical basis is that 
when the basis vector changes direction, it can pick 
up components along only the e^-i and e^+i directions. 
Furthermore, because ^(e^+i • e^) = 0, we have the 
relation [23l 



dN dN 



Hence, the matrix 



Den 

dN 



(24) 



(25) 



is skew-symmetric with the only nonzero components be- 
ing 



An+1 
(n) ' 
Xn 



(26) 



The Z matrix therefore encodes how quickly the d basis 
vectors are changing direction. For comparison, we note 
that the derivatives of the basis vectors and the Z matrix 
are equivalent to those in [23 , even though we defined 
the kinematical vectors differently. 

Lastly, we compare Z when there are one, two, and 
three or more fields present during inflation. In single- 
field inflation, there is only one basis vector and this basis 
vector cannot change direction, so Z reduces to the zero 
matrix. In two-field inflation, the only nonzero coeffi- 
cients of Z are Z21 = —^12 = ^5 and so only this single 
quantity is needed to characterize how quickly both basis 
vectors are covariantly changing directionj^ In [13 , we 
called this quantity the turn rate^ since it represents how 
quickly the background field trajectory is turning with 
respect to the field manifold. However, for multi-field in- 
flation with three or more fields, there are multiple such 
rates that characterize how quickly the d basis vectors are 
rotating. Therefore, to extend our two-field formalism to 
multi-field infiation, we call | \ the turn rate for the 
nth basis vector^ and we say that the kinematical quan- 
tity Zn+i,n represents how quickly the basis vector is 
turning into the direction of e^+i. Because Z summa- 
rizes the turn rates for all basis vectors, we call Z the 
turn rate matrix. The turn rate matrix is therefore the 
multi-field generalization of the idea of a single covariant 
turn rate for two-field infiation. The turn rate matrix, 
along with the kinematical basis vectors, plays a key role 
in determining the dynamics of the field perturbations. 



C. Slow-Roll And Other Limits 

Having delineated a framework to represent the kine- 
matics of any multi-field model, we now consider the 
slow-roll limit and other special subcases of limiting be- 
havior that simplify the equations of motion. 

In standard single-field infiation, the slow-roll limit is 
invoked in order to guarantee nearly exponential expan- 
sion of the scale factor, which produces a nearly scale- 
invariant density power spectrum. In this limit, the field 
is slowly rolling, which effectively means that both the 
field speed and field acceleration are small, satisfying 
77 <C <C 1. (Recall that we constructed both the field 
velocity and field acceleration to be dimensionless.) The 
slow-roll limit implies that the field acceleration is negli- 
gible in Eq. ^ , which allows us to simplify the equations 
of motion and obtain analytic estimates for the power 
spectrum and spectral observables. 

In multi-field infiation, an analogous slow-roll limit is 
typically defined by the two conditions (e.g., [29 } [32| \33 \ 



^ In [13], we wrote the quantity ^ as ^ since in two-field infia- 
tion, there is only a single direction perpendicular to the field 
trajectory. 
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1351 EH): 



e«i|VlnT/|'«l, 



|M,,-|«1. 



(27) 
(28) 



In some approaches (e.g., [TH [HJ [23] ) , the second condi- 
tion above is effectively replaced by the condition 



(29) 



in our notation. We examined these conditions in the 
case of two-field inflation, the simplest case of multi-field 
inflation, in [T3^. As a result of that analysis, we argued 
for a more nuanced approach to the slow-roll approxi- 
mation that splits the slow-roll limit into two separate 
limits — the slow-roll limit and the slow-turn limit. Our 
redefined slow-roll limit corresponds to the conditions 



e < 1, 
«1, 



(30) 



which means that the field speed and the rate at which 
the field speed is logarithmically changing are small. This 
redefined slow-roll limit corresponds to limits on single- 
field-like behavior. In contrast, our two-field slow-turn 
limit in [13] corresponds to what can be viewed as dis- 
tinctly multi-field behavior. It corresponds to limits on 
how quickly the field trajectory is covariantly changing 
direction — a distinctly multi-field behavior — and it can 
be expressed as 



Dei 



dN 



Z21 <c 1. 



(31) 



The power of distinguishing between the rolling and turn- 
ing behaviors of the field vector is that they have very dif- 
ferent effects on the density power spectrum, and differ- 
ent approximations can be made depending on whether 
the field velocity is slowly rolling, slowly turning, or both. 

Here, we extend this more nuanced approach for two- 
field inflation to general multi-field inflation. First, we 
will refer to the redefined slow-roll limit — the limits on 
single-field-like behavior, as defined by Eq. (30) — as the 
single-field slow-roll limit Second, to extend the slow- 
turn formalism to multi-field inflation, we say that a basis 
vector is slowly turning if 



Den 



dN 



«1, 



(32) 



and also say that the corresponding scenario is in the 

slow-turn limit. By Eq. (bsL if Xn-i^ •> Xn \ ^"^^ 

Xn-^i^ nonzero, then Eq^32| is equivalent to the 

requirement that X^n-i^ ^ Xn^^ ^ X^n-^i^ • When all d 
basis vectors are slowly turning, we say that the inflation- 
ary scenario is in the slow-turn limits and the magnitude 
of every component of the turn rate matrix, Z, is sig- 
nificantly less than one. For reasons that will become 



more apparent later, the full slow-turn limit represents 
limits on multi-field behavior; this is because the turn 
rates play a critical role in determining the interactions 
among modes in the kinematical basis. 

Taken together, the single-field slow-roll limit and the 
slow-turn limit constitute our two main kinematical lim- 
its of interest. However, in investigating the modes in 
the kinematical basis, we will want to invoke not just 
these two limits but also the limit in Eq. (28) to en- 



sure that the acceleration term in the mode equations 
can be neglected. When we combine these three limits 
together, we call the combined limit simply the slow-roll 
limits even though we are also assuming full slow-turn 
behavior. We do this to avoid introducing new nomen- 
clature that would create confusion with past literature. 
In this paper, we focus on the full slow-roll limit, but 
we introduce the above separate limits so that different 
individual limits can be explored, like we did for two- 
field inflation in [13]. For example, by exploring these 
separate limits in [13 , we found that two-field models 
that strongly violate the slow-turn limit around horizon- 
crossing but not the single-field slow roll limit are ruled 
out by WMAP constraints on the power spectrum. We 
emphasize that much is to gain from a similar thorough 
investigation of multi-field inflation, and this paper lays 
much of the foundation for such an investigation. 

Now we are ready to apply the above set of limits to 
the background equations of motion and our kinematical 
vectors to simplify them. Eq. ([t]) for the evolution of the 
fields reduces to 



(33) 



whenever the single-field slow-roll and the ei slow-turn 
limits hold; the full slow-roll limit does not need to be 
invoked. Rotating to the kinematical basis, in the single- 
field slow-roll and ei slow-turn limits, the field speed is 
given by 



or equivalent ly by 



iVlnVl 



(34) 



(35) 



Also under the same two limits, by virtue of Eq. (33), 
the operator ^ = 0' ' V becomes 



D 

dN 



-VlnF- V. 



(36) 



Using this expression, the kinematical vectors in the full 
slow-roll limit can be approximated by 



An) 



For example. 



{-VlnV- Vy^-^^ (-V^ InV). (37) 



T] ^ -Met)' ^ MVUnV, 



(38) 
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and 



-Mr] 



DM 

dN 



0' 



[M^ + ( V In y • VM)] In V. 



(39) 



The slow-roll approximations for the kinematical basis 
vectors follow directly from the above results, with Eq. 
(37) substituted for x^^^ in Eq. (16). For example, in 



the slow-roll approximation. 



ei 



e2 



VUnV 
|VlnV|' 

V^lnV 



VlnV 



|VlnV| V|VlnV| 



M In V, 



(40) 



and so on. 

These results will be very important in simplifying the 
equations of motion and the interactions among modes 
in multi-field inflation. 



III. FIELD PERTURBATIONS 

Now we consider the field perturbations during in- 
flation. First, we present equations of motion for the 



field perturbations in Section |III A[ We provide evolu- 
tion equations in both the given and kinematical bases. 
In Section |IIIB[ we explore the properties of the mass 
matrix, showing for the first time that most of its coef- 
ficients vanish in the kinematical basis. We use this re- 
sult in Section [ill C| to show that the interactions among 
modes greatly simplify in the kinematical basis. We then 
analyze these mode interactions and discuss how they de- 
pend on the geometric features of multi-field inflationary 
Lagrangians and their associated field trajectories. This 
section lays the key foundation that allows us to connect 
the geometric features of inflationary Lagrangians to the 
cosmic observables. 



A. Evolution Equations for the Field Perturbations 

In this section, we present equations of motion for the 
field perturbations. In doing so, we work in the flat 
gauge, and we do this for two reasons. First, in this 
gauge, the field perturbations decouple from the met- 
ric perturbations. Second, the field vector perturbation 
in the flat gauge turns out to be the gauge- invariant 
Mukhanov- Sasaki variable, 6<pf = Scj) -\- tljcj)\ where 
represents the scalar metric perturbation on spatial hy- 
persurfaces [63, 64 . From here forward, since we con- 
sider the perturbations only in the flat gauge, we drop 
the subscript / from Scj)f. 

We derive an equation for the field perturbations in the 
given basis by imposing covariant conservation of energy 
to first order in the field perturbations, as has been done 



before [TTl [181 [21] . However, instead of the standard ap- 
proach, we use as the time variable, both because it 
is more physically intuitive and it makes the equation of 
motion dimensionless. We derived such an equation in 
the context of general two-field inflation with noncanon- 
ical kinetic terms in [13]. Following the same series of 
steps, we arrive at the same expression, with the excep- 
tion that the curvature term arising from the field metric 
is a matrix in general multi-field inflationj^ The resulting 
equation in Fourier space is 



(3- 



1 D^Scj) DScj) 

dN 



dN^ 



M 



64) 



(3-e)2 



50, (41) 



where k is the comoving wave number. The term M is 
the effective mass ma^raj^and we define it as 



M = M 



1 



(3-6) 

and the curvature matrix, R, is defined as [21] 



(42) 



(43) 



where i^^^^d Riemann curvature tensor associated 

with the field metric. Because of the symmetry and anti- 
symmetry properties of the Riemann curvature tensor, it 
follows that R is symmetric. Moreover, R(/)' = 0, and 
hence R = R^x. 



Next, we consider various limits of Eq. (41). In the 



superhorizon limit, when the modes are significantly out- 
side the horizon such that (^) <C 1, the subhorizon 
term ^^ly ^ Scj) can be neglected. In the combined single- 
field slow-roll and ei slow-turn limits, we argued in [13] 
that the term (^z^ can be neglected in two-field infla- 
tion since this term is much smaller than IVI. The exact 
same arguments as in [13 hold for general multi-field 
inflation with an arbitrary number of fields. And, when- 
ever both the superhorizon limit and the full slow-roll 
limit apply, it can be shown [23l [65] that the acceleration 
of the field perturbations can also be neglected. In these 
combined limits, Eq. (41) reduces tc[^ 



D5(j) 

~dN' 



-MScf). 



(44) 



^ For comparison, in two-field inflation, the curvature term effec- 
tively reduces to a single degree of freedom: the Ricci scalar 
(times a scaled outer product of two kinematical basis vectors). 

^ In comparison to our dimensionless deflnition of the effective 
mass matrix, Nibbelink and Van Tent deflned the effective mass 
matrix as = VV^y - H'^H [23]. 

5 Technically, Eq. ^ also assumes that the dimensionless com- 
ponents of R are much less than one. 
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Let us consider the above superhorizon, slow-roll equa- 
tion of motion in the kinematical basis. In the kinemat- 
ical basis, there are up to d nonzero field perturbations, 
and the nth mode is 



en • Scf). 



(45) 



The benefit of working in the kinematical basis is that it 
allows us to separate out the adiabatic mode from the en- 
tropy modes. The adiabatic or density mode corresponds 
to S(j)i , meaning that it corresponds to the component of 
Scj) that is parallel to the field trajectory. The remaining 
components of ^0 in the kinematical basis correspond to 
entropy modes. Entropy modes are linear combinations 
of the field perturbations that leave the overall density 
unperturbed. As there are d fields in the system, there 
will he d — 1 entropy modes, all of which are orthogonal 
to the field trajectory and to each other. 

Now the equation of motion assumes a different form 
in the kinematical basis because the basis vectors can 
rotate, causing the equation of motion (41) to pick up 



extra terms that trivially vanish in the original given ba- 
sis. Using the fact that 



(46) 



it follows that the corresponding equation in the kine- 
matical basis is 



DS(t) 

"dW 



M + Z 



(47) 



Note that the term in brackets is the sum of a symmetric 
matrix and an anti-symmetric matrix. If, however, we 
are interested in just the evolution of the d — 1 entropy 
modes in the kinematical basis, their equation of motion 
in the combined superhorizon and slow-roll limits is 



dN 



M 



Set) 



±5 



(48) 



where the special matrix projection was defined in Eq. 
(20). Therefore, if we want to determine the evolution of 



only the entropy modes, we can study them by reducing 
the dimension of the perturbed equation of motion by 
one. 



B. Properties of the Mass Matrix 



([t]). Projecting this equation onto the basis vectors gives 
• V In V = for n > 2, since in the kinematical basis, 
T] has nonzero components along only ei and 62- Con- 
tinuing, we differentiate the equation of motion for the 
fields to obtain Eq. (12) for the jerk, ^. Projecting this 
equation onto the basis vectors gives Min = for n > 3, 
since ^ has nonzero components along only ei, 62, and 
63. Similarly, repeatedly differentiating Eq. ^ and pro- 
jecting the result onto the basis vectors allows us to find 
the coefficients of M and its covariant derivatives that 
vanish. In the end, we find that 



Mmn for \m-n\> 2, 



for |m-n| >p + 2. (49) 



The second condition above implies that the vector 
formed from contracting V^V^ with n — 1 instances of 
(p^ vanishes along the em direction, where m > n + 1. 
This shows that most of the coefficients of M vanish in 
the kinematical basis. 

In the full slow-roll limit, the mass matrix simplifies 
even further. The proof is similar to the above proof 
for the general case, except we start with the slow-roll 
equation of motion ( [33| ). Projecting this equation onto 
the basis vectors trivially gives that • V In V = for 
n > 1. Differentiating the slow-roll equation of motion 
once and projecting it onto the basis vectors gives Min = 
for n > 2. Similarly, repeatedly differentiating Eq. 
(33), we obtain 



Ari+2) 



■E 

m=0 



D 

dN 



M 



^(n— m+l) 



• (50) 



Projecting Eq. (50) onto the kinematical basis vectors 
gives 



D 

dN 



M 



for \m — n\ > 1, 

for |m — n| > p + 1. (51) 



Therefore, in the full slow-roll limit, M reduces even fur- 
ther and becomes a symmetric tridiagonal matrix in the 
kinematical basis. This concludes our proofs showing 
that most of the coefficients of the mass matrix do in- 
deed vanish in the kinematical basis. 



In general, the interactions among modes are most eas- 
ily understood in the kinematical basis. First, in this ba- 
sis, the density mode can be teased out from the d — 1 
entropy modes. Second, it turns out that most of the 
coefficients of the mass matrix vanish in this basis. This 
second fact profoundly simplifies the interactions among 
modes, as we will demonstrate in the following section. 

But first, in this section, we prove our claim that most 
of the coefficients of M vanish in the kinematical basis. 
We start from the full equation of motion for the fields 



C. Mode Sourcing Equations 



Our goals in the rest of Section (III) are two-fold: to 
simplify the equations of motion as much as possible and 
to gain intuition and physical insight into the mode in- 
teractions. Doing so will bring us much closer to our 
ultimate goal of connecting the geometric and kinemati- 
cal features arising from the inflationary Lagrangian with 
the cosmic observables. 
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In exploring the interactions among modes, we con- 
sider the superhorizon equation of motion for each of the 
d modes in the kinematical basis. When a mode 5(j)m af- 
fects the evolution of mode S(j)ri^ we say that 5(j)rn sources 
S(j)n, regardless of whether that interaction causes d(j)n to 
increase or decrease in amplitude. And we call the in- 
teraction between the two modes a sourcing relationship. 
Since we will group all sourcing terms on the right-hand 
side of each equation, we refer to such equations as mode 
sourcing equations. Surprisingly, studying the interac- 
tions among modes one by one like this has not received 
much attention in general multi-field infiation; however, 
studies have been undertaken concerning the adiabatic 
mode in multi-field infiation (in particular, [2TJ[23]) and 
in general two-field infiation [131 129H4Q] . Here we fill this 
important gap in the literature. 

We start with the mode sourcing equation for the adi- 
abatic mode, 6(f)i^ which is most easily derived from the 
fact that the comoving density perturbation vanishes in 
the superhorizon limit. Imposing this constraint yields 

[El ESI [30] 



^ =2^.^ = 2Z.,(;^). (52) 
V J dN V \ V ' 



In the slow-roll limit, the above equation can be written 
as 



•Mil 



2Z, 



2lC 



(53) 



where we have used that Z21 ~ — M12, which follows 
from Eqs. (I2TI), (l38|), and (l52|. Examining Eq. §^ 



the field trajectory changes direction with respect to the 
field manifold (e.g., [1 [H [H] [H [30]). The strength 
of the sourcing depends on the ei turn rate, so in our 
kinematical picture, the faster the background trajectory 
changes direction, the more the adiabatic density mode 
grows Moreover, the adiabatic mode can be sourced 
only by the entropy mode S(j)2'-, none of the other modes 
can source the adiabatic mode. Otherwise, when the field 
trajectory does not turn or S(j)2 vanishes, it follows that 
6(j)i cx V, which is tantamount to single-field behavior. 

Continuing, for each of the entropy modes, we can also 
derive a mode sourcing equation. Here we present the 
results only for the slow-roll limit, as the full equations 
of motion are significantly more cumbersome without of- 
fering much additional insight, and moreover, they are 
usually not needed. For the d—1 entropy modes, we can 
derive slow-roll approximations to the sourcing equations 
starting from Eq. (47) and using Eq. (51). For the 6(j)2 
mode, we obtain 



22C 



1 

(M23 - ^32) Hs^^Yl ^2m^ 



(54) 



reveals that the adiabatic mode is sourced only when 



where we have again used that Z21 ~ — M12. First, notice 
that the adiabatic mode does not source the d(j)2 mode; 
in fact, the adiabatic mode does not source any of the 
entropy modes in the superhorizon limit (regardless of 
whether the slow-roll limit applies). Second, if the kinetic 
terms are canonical, then R vanishes, and (5^2 can be 
sourced only by S(j)s. Conversely, S(j)2 will be sourced by 
the S(j)n mode whenever i?2n 7^ 0- Similarly, for the S<pn 
mode, where n > 3, the sourcing equation is 



,n-l + ^n,n-l j ^^n-l - [Mn 

m=2,|n— m|>l 



Importantly, the coefiicients of the mass and turn rate 
matrices allow 5(j)ri to be sourced by only two modes: 
3(j)n-i and (^(/)n+ij^ If the kinetic terms are noncanoni- 
cal, the other entropy modes 5(j) Yn can source (50^ only if 



In the classical treatment of Eq. ([52|, this statement is implied 
to be taken with respect to assuming that and 5(f)2 are both 
positive. In the quantum treatment, when we speak of the mode 
5(1)1 growing, it is implied that we are referring to the variance 
of (5^1 growing. Similar assumptions are made when discussing 
the other modes. 

For comparison, in the non-slow-roll case (results not shown), 
up to four modes ((5(/)n-2, ^(/'n-i , ^</>n+i , and (5(/)n+2) can source 
5(t)n when the kinetic terms are canonical. This follows directly 
from the special form of the mass matrix in the kinematical basis, 
which we covered in the previous section. 



Rnm 7^ is satisfied. These results are tantamount to re- 
ducing the complexity of the matrix equation (47) by rec- 



ognizing that the sum M+Z is a tridiagonal matrix in the 
kinematical basis and by recognizing that M12 = —^21- 
These novel results show that the sourcing equations 
for the modes profoundly simplify in the kinematical ba- 
sis, which will allow us to better understand and more 
easily calculate the interactions among modes. 



1. Connections between the mode sourcing and the 
geometry and kinematics of inflation 

Now we analyze the connections between the mode 
sourcing and the geometry and kinematics of infiation. 
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First, consider what happens in the absence of sourc- 
ing. For 5(j)n to be unsourced, Eqs. (53)-(55) show that 



must not be turning and the mass matrix coefficients 
^n,n+i and Mn,n-i must vanish, as well as the curvature 
matrix coefficients Rnm foi" all m 7^ n (or when the sourc- 
ing terms in aggregate sum to zero). When these condi- 
tions are met, the 5(j)ri mode is unsourced and obey^ 



- f Mr^ridN 



(56) 



That is, the effective mass for the nth mode. Mnn^ deter- 
mines that mode's intrinsic evolution rateO 

The effective mass depends on the covariant Hessian 
of the inflationary potential, Mnn = VnVn In and on 
the curvature tensor of the field manifold. Since these are 
both geometric quantities, this means that we can predict 
the behavior of (5^^ simply by determining the geometric 
features of the Lagrangian. For instance, when the po- 
tential In V is "concave up" along the direction — that 
is, when the Hessian coefficient Vt^Vt^, InV is positive — 
this causes 6(j)n to decay. Conversely, if the potential 
has a concave down shape along the direction, this 
causes to grow. A well-known example exhibiting 
this behavior is the adiabatic mode in single-field infla- 
tion, which is often likened to a ball rolling down a hill 
that speeds up or slows down depending on the ground's 
concavity. A second good example is the 5(j)2 mode in 
two-field inflation with canonical kinetic terms; if the 
two-dimensional field trajectory is rolling in a valley (con- 
cave up), then the S(j)2 mode decays, but if it is rolling 
along a hill (concave down), then 5(j)2 grows in amplitude. 

The second geometrical quantity, the curvature term 
Rnn^ is a bit more complicated to understand. It involves 
the contraction of the Riemann tensor of the field metric 
with two instances each of ei and e^. Geometrically, it 
represents 2e times the component of the failure of ei 
to be parallel-transported around a closed loop defined 
by the directions ei and e^. Therefore, its contribution 
to the mode sourcing can in principle be determined from 
the geometry of the field metric and the potential. If this 
deviation from parallel transport of ei results in a posi- 
tive component along the direction, then this causes 
6(j)n to grow; conversely, a negative value causes S(j)n to 
decay. For example, in two-field inflation, since R22 is 
proportional to negative e times the Ricci scalar of the 
field manifold, if the field manifold is locally elliptical, 
this will cause S(j)2 to decay, while a locally hyperbolic 



^ We caution the reader that the mode amphtude should not be 
assumed to decay approximately as e~^rin^* ^ where M*^ is the 
value of the effective mass at horizon exit and A/"* is the number 
of e-folds since the mode exited the horizon. As we discussed in 
[13| for the case of two-field infiation, this assumption often leads 
to large inaccuracies in estimating the mode amplitudes and the 
spectra. See 13 and references therein for further discussion of 
this issue. 

^ However, we note that very often the curvature term is smaller 
than the mass M^n, and so frequently, Mnn ~ Mnn- 



surface will cause 6(j)2 to grow (see fT3^). Note that if 
both Mnn and Rnn have the same sign, they will par- 
tially negate each other; thus, it is their net effect that 
matters. Therefore, we can view the effective mass as 
some sort of measure of the "net curvature" or geometry 
of the inflationary Lagrangian along the single direction 
specified by e^. And the steeper this net curvature of the 
Lagrangian along the direction, the faster the mode 
S(l)n evolves. 

The corollary of Eq. (56) is that when S(j)n is un- 
sourced, the quantity 



(57) 



is conserved in the superhorizon limit. For example, in 
single-field inflation, the S<j)i mode is automatically un- 
sourced, and hence the quantity 



V 



(58) 



is conserved. In inflation with two effective fields, the 
entropy mode 6(^2 is unsourced, so the quantity 



t>2e 



J M22dN 



(59) 



is conserved. (Thus, in the special case of two-field in- 
flation, the above expression for S(j)2 can be plugged into 
Eq. ( 53 ) to find a semianalytic expression for , elimi- 



nating the need to solve a set of coupled equations [13].) 
Therefore, we make the important observation that Eq. 
(5^ is the multi-field generalization of the well-known 
single-field conservation equation for the adiabatic mode 
in Eq. (58). Eq. (57) endows each of the d modes with 



a conservation equation that holds whenever is not 
turning, M^^^-i Mn^n+i 0, and Rnm for ah 
m ^ n (or when the sourcing terms in aggregate sum to 
zero). So there are up to d potential conserved quantities 
related to the modes. 

Next, we consider how sourcing affects the evolution 
of the 5(j)n mode. According to Eq. (55), sourcing effects 
can arise from: 

1. Off-diagonal terms in the mass matrix 

2. Any nontrivial geometry of the field manifold 

3. The kinematical basis vectors changing direction 

We will consider each sourcing effect in turn. 

The first type of sourcing effect is a geometric ef- 
fect that arises from off-diagonal terms in the covari- 
ant Hessian, M, of the inflationary potential. Specifi- 
cally, the mass matrix M introduces terms of the form 
— Mn,n-i^^n-i and — Mn+i,n^^n+i mto the sourcing 
equation for d(j)n' These terms can be viewed as mea- 
sures of the coupling between fields in the potential and 
of whether this coupling results in a higher or lower po- 
tential energy state. Or in the geometric picture, the 
terms Mn^n-i and Mn^n+i represent how much the nth 
component of the covariant derivative of In V varies along 
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the Bn-i (Mn,n-i) and e^+i (M^^n+i) directions. If 
the couphng term Mn,n+i (or Mn,n-i) is positive, then 
S(j)n-\-i (^^(/>n-i) win cause S<pn to decay; otherwise, if it 
is negative, it wih increase the amphtude of S(j)ri- And 
the greater these geometric features in the potential, the 
stronger the sourcing effects will be. Interestingly, since 
the mass matrix is symmetric, a nonzero Mn^n-\-i leads 
to parallel sourcing effects; for example, a negative value 
for Mn^n±i will cause both the 5(j)ri and 5(j)ri+i modes 
to growj^ Therefore, knowing either the mathemati- 
cal form of the inflationary potential or the potential's 
geometric shape provides insight into the nature of the 
sourcing effects. 

Likewise, the geometry of the field manifold can also 
cause sourcing effects. As explained earlier, the form of 
the kinetic terms in the inflationary Lagrangian can be 
represented through a field metric, and this metric can 
be viewed as inducing a field manifold. If the field man- 
ifold has nontrivial geometry, then the Riemann curva- 
ture tensor will be nonzero, and this will be manifested 
in the form of a nonzero symmetric curvature matrix 
= 2eRl^^e\e\. Specifically, if the component of 
the failure of ei to be parallel transported around the 
closed loop defined by ei and is nonzero, then the 
curvature matrix will cause 5(j)rn to source 5(j)n' Since 
the curvature matrix can be factored into e times a term 
involving the Riemann tensor, this term technically com- 
bines geometric and kinematical effects; so when all else 
is equal, the impact of noncanonical terms on the mode 
sourcing tends to be greatest at the end of inflation and 
whenever else the field speed is large. For this reason, 
the curvature matrix may be better viewed as a measure 
of the impact of nontrivial kinetic terms on the mode 
sourcing. However, for simplicity, we tend to refer to 
this type of sourcing effect as a geometrical effect arising 
from the field manifold. Now like the mass matrix, since 
the curvature matrix is symmetric, a positive value for 
a given curvature matrix coefficient will cause both the 
5(j)rn and 5(j)n modes to grow. However, in comparison to 
the sourcing effects due to the mass matrix coefficients, 
there are two key differences: (1) the curvature matrix 
can in principle couple together any two modes, and (2) 
the curvature matrix appears in the equation of motion 
with the opposite sign. Thus, we may view the mass ma- 
trix and curvature matrix as together encapsulating the 
sourcing effects due to the geometry of the Lagrangian, 
with the mass matrix primarily corresponding to the po- 
tential and the curvature matrix to the kinetic terms. 

The third and last kind of sourcing effect is purely a 
kinematical effect — a direct consequence of the kinemat- 
ical basis vectors changing direction. Consider the term 
Zn+i,n^(t^n+i in Eq. ( [55| ). The kinematical term Zn+i,n 
represents how quickly the basis vector is turning into 



Sourcing Terms 


What the Terms Represent 




Covariant Hessian of Potential 
(geometry of potential) 


Rnm 


T ■) • rn £ T^' 1Jl\/r T1J 

Riemann lensor or lieid Maniroid, e 
(geometry of field manifold, kinematics) 




Turn Rate of (kinematics) 



TABLE I. The three types of sourcing effects in the mode 
sourcing equation for 5(j)n and what each set of terms effec- 
tively represents. More detailed explanation of the terms and 
their impact on mode sourcing is given in the text below. 



the direction of the e^+i basis vector. Since Z^+i^^ is 
always non-negative, this turning will always cause 5(j)n 
to grow. And the faster is turning into the e^+i direc- 
tion, the more 5(j)n+i sources 5(j)n' This sourcing effect 
can be interpreted physically as follows: the direct rota- 
tion of the kinematical basis vectors causes what was once 
a 5(j)n+i mode to be partially converted into a 5(j)n mode. 
The other kinematical sourcing term, —Zn^n-i^^t^n-ii can 
be understood similarly. However, this term causes 5(j)ri 
to shrink in magnitude, which can be explained by the 
fact that 5(j)n is being partially converted into by 
the rotation of basis vectors. The anti-symmetry of the 
turn rate matrix neatly encapsulates these antithetical 
kinematic effects. 



A natural question that arises is whether the kinemat- 
ical terms in the turn rate matrix represent some geo- 
metric feature of the Lagrangian. The answer is that 
while they can be expressed in terms of V and the field 
metric via Eq. (37), this expression is not geometrically 



Again, when working in the classical picture, our statements are 
with respect to positive field fiuctuations; it is straightforward 
to extrapolate to other cases. 



transparent and does not boil down to a simple geomet- 
ric feature. For example, one might expect that the mass 
matrix somehow effectively determines all the turn rates, 
but in general it does not, and the turn rates involve 
more complex combinations of the various nth-order co- 
variant derivatives of In y . The only common exception 
occurs in the case of quadratic potentials with ca nonical 
kinetic terms (which is covered in Section HI C 3[ ), where 
it turns out that all the turn rates can indeed be ex- 
pressed in terms of coefficients of the mass matrix. This 
leads to a very special simplification in the mode sourc- 
ing equations and also offers new insight into multi-field 
inflation. But in other scenarios, it is often easiest to 
view the effects from the turn rate matrix as kinematical 
effects, rather than particular geometric effects. 
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2. Special Case I: Canonical kinetic terms 

When the kinetic terms are canonical (or less strin- 
gently, when R is no more complex than tridiagonal in 
the kinematical basis), the sourcing equations simplify 
further. In this common case, the mode sourcing equa- 
tion for the adiabatic mode is unchanged in form. How- 
ever, the slow-roll mode sourcing equations for the en- 
tropy modes simplify to 



22<~ 



(M23 - ^32) ' 



(60) 



and 



- (^Mn,n+1 - ^n+l,n) (61) 



for n > 3. Therefore, when the kinetic terms are canoni- 
cal, only up to two modes, 5(j)n-i and ^^n+i, can source 
the 3(j)n mode in the slow-roll limit. Solving this series of 
coupled equations for the modes is tantamount to solv- 
ing slow-roll Eq. (47) where the matrix sum M + Z is 
tridiagonal. 



Special Case II: Quadratic potentials with canonical 
kinetic terms 



When the kinetic terms are canonical and the potential 
is quadratic — meaning that it satisfies 



WW = 



(62) 



— then the mode equations simplify even further. Since 
= for n > 3, it follows that (^)^M, where 
p > 0, is a function of only VlnV and M. Therefore, 
projecting the kinematical vectors in Eq. ([50| onto the 
basis vectors yields 



(n+1) 
X(n) 



X 



(n) 



n+l,n • 



(63) 



That is, the mass matrix coefficient — M^+i^^ gives the 
rate at which the basis vector turns into the direction 
of the e^+i basis vector. Substituting this result into Eq. 



(61), the mode sourcing equation for all modes reduces 



to 



(64) 



Therefore, for quadratic potentials with canonical kinetic 
terms, the d(j)n mode is sourced only by the mode 
and only when the basis vector is covariantly changing 
direction into the e^+i direction. Mathematically, the 
solution for the nth mode is 



N 



2Z, 



(65) 



where * denotes a quantity is to be evaluated at horizon 
exit. In principle, one could solve the above series 
of integrals. Finally, since there is no sourcing when 
^n+i,n = the number of kinematical basis vectors that 
are changing direction inversely indicates the number of 
conserved mode quantities. 



4. Prescription for analyzing how Lagrangian differences 
affect mode sourcing 



Taking a closer look at the special result in Eq. (63) 



for quadratic potentials, this result follows from the more 
general result in Eq. ( 50 ) that 



^n+l,n ^ -Mn+l,n + COrrCCtioUS. 



(66) 



These so-called "corrections" depend on certain projec- 
tions of the higher-order covariant derivatives of \nV] 
that is, each term contains at least one quantity V^F, 
where p > 3. The highest-order term in the series of cor- 
rections involves V^^F, and the sum of correction terms 
truncates prematurely if V^F vanishes for all p greater 
than some positive integer but less than n. Interestingly, 
plugging Eq. (66) into the general mode sourcing equa- 
tion (55) tells us that the sourcing of by S(j)n-i is 
controlled by these certain projections of the higher-order 
covariant derivatives In where p > 3. Similarly, we 
can view the sourcing of 5(j)ri by S(j)n-\-i to be controlled by 
a term that is twice the turn rate Z^+i^^, plus corrections 
from higher-order covariant derivatives of InV. And of 
course, any nonzero Rmn terms arising from noncanon- 
ical kinetic terms would contribute additional sourcing 
effects. In other words, we can view the entropy mode 
sourcing equations as 



/)^ + MnnS(l)n ^ — (corrcctious from potential) (^^n-i + (2^n+i,n + corrections from potential) ^^n+i 
1 ^ 

(67) 
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where the "corrections" are certain projections of the co- 
variant higher-order derivatives of the potential, as ex- 
plained above. This results in a very interesting and 
useful way to view the interactions among modes: the 
interactions can essentially be divided into sourcing ef- 
fects shared in common with canonical quadratic models 
(Zn-\-i,n terms) and sourcing effects arising from devi- 
ations from this fundamental Lagrangian (the ^'correc- 
tions^' and Rnm terms). We advocate this novel approach 
as a powerful prescription for exploring how differences in 
Lagrangians translate into differences in mode dynamics. 

This completes our more rigorous and physically- 
motivated study of the interactions among modes in gen- 
eral multi-field inflation. 



number of these basis vectors gives the dimension of the 
background solution. Because of the way we constructed 
the kinematical basis vectors in Eq. (16), the dimension 



IV. 



EFFECTIVE NUMBER OF FIELDS AND 
THE SPECTRAL OBSERVABLES 



In Sections [II| and IIL we explored how the inflationary 
geometry and kinematics affect the evolution of modes. 
In this section, we discuss how this in turn translates into 
the cosmic observables. But since determining how the 
geometry and kinematics affects the mode interactions is 
most of the battle, we will not rehash any connections 
that can easily be inferred from our discussion in Section 
III CI Instead, here we will focus more on how the infla- 



tionary geometry and kinematics determine the effective 
number of inflationary fields and how this is reflected in 
the cosmic power spectra, bispectrum, and trispectrum. 
We start in Section |IV A| by discussing how the kine- 
matics and geometry of inflation determine the effective 
number of fields. In Sections IV BpV Pj we discuss the 
multi-field results for the power spectra, bispectrum, and 
trispectrum. In the process of doing so, we introduce a 
new cosmic multi-field observable that distinguishes two- 
field models from models with three or more fields, and 
we present a new multi-field consistency relation. 



A. Effective Number of Fields 

Having laid the groundwork for determining the sourc- 
ing relationships among the modes, we can now consider 
the effective number of fields in a mult i- field scenario. 
We define the effective number of fields to be the mini- 
mum number of fields necessary to adequately describe 
an inflationary potential and reproduce the results for the 
cosmic observables across the distance scales of interest. 
The effective number of fields can also be thought of as 
the dimension of inflation. 

We start by considering the minimum number of fields 
required to produce a given background solution. In or- 
der to adequately describe the background solution, we 
must have enough fields to reproduce all the kinemati- 
cal vectors, as defined in Eq. (14). That is, we must 



also equals one plus the number of kinematical basis vec- 
tors that are changing direction. In particular, if the ei 
basis vector is changing direction, then the inflationary 
scenario has at least two effective field degrees of freedom. 
Of course, at a given moment, the effective field degrees 
of freedom may be fewer or greater than it is at other 
times. To determine an effective dimension for an entire 
multi-field scenario, one can simply take the maximum 
dimension of the kinematical vectors that exists between 
the start of that scenario until the end of inflation. This 
is illustrated in Figure [l] In scenarios with canonical ki- 
netic terms that are effectively described by one field, the 
trajectory in field space will resemble a line. Similarly, 
the trajectories of scenarios with two effective fields will 
reside in a plane. 

The effective dimension of the field perturbations, how- 
ever, is more complicated to determine. We define the 
effective dimension of the modes as the minimum num- 
ber of independent field perturbations we must consider 
in order to calculate the spectra. Based on the spectra 
we consider, this is equivalent to the number of fields 
we must consider in order to find solutions for 5(j)i and 
5(j)2] why this is so will become clearer in the next two 
sections. 

Consider first the particular case where the ei basis 
vector never turns. In this case, the adiabatic mode is 
never sourced and the quantity ^ is conserved in the 
superhorizon limit. While the effective dimension of the 
unperturbed fields is one, the dimension of the perturbed 
fields can be either one or more: it is one if there is only 
one field and hence there are no entropy modes, but it is 
more if there are two or more fields during inflation and 
hence a power spectrum of entropy modes. 

In the case where the field trajectory does change di- 
rection during inflation, there are two reasons why the ef- 
fective dimension of the background and perturbed fields 
do not necessarily coincide. The first reason is that the 
curvature matrix R can couple together the various en- 
tropy modes, independently of the turning behavior of 
the kinematical basis vectors. Second, in general, it is not 
true that M^+i^^ ^ —Zn^n+i- So even if the kinematical 
basis vectors are not turning, a nonzero Mn+i,n could 
still allow the ^(/)n+i mode to source the mode. (Of 
course, for many models, when Z^+i^^ = 0, it will also 
be true that Mn+i^n = 0-) Therefore, for models with 
at least two fields, the effective number of field pertur- 
bations we need to consider in order to find expressions 
for S(j)i and S(j)2 will equal two plus the number of con- 
secutive nonzero terms Mn+i,n + ^n+i,n when we count 
upwards starting from n = 2. This follows directly from 
the series of slow-roll sourcing equations in Eq. (55). 



have enough kinematical basis vectors to span the space 
defined by the kinematical vectors, and the minimum 



Based on this analysis, we take the effective dimension 
of the perturbed field system, which can be larger than 
the dimension of the unperturbed system, as the overall 
effective dimension of a multi-field scenario. 
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Yet, although we can assign an overall dimension to 
each scenario, it is also useful to consider that an in- 
flationary scenario may be broken into multiple phases, 
with each one defined by a different effective number of 
fields being active. For example, in canonical quadratic 
models with very different masses for the fields in the po- 
tential, there are periods dominated by the dynamics of 
a single-field, punctuated by periods in which two fields 
dominate the dynamics. By understanding that a model 
with multiple fields can be approximated by a series of 
scenarios with a much smaller effective dimension — in the 
example, a series of single-field and two-field scenarios — 
we can gain much greater insight into the key features of 
such models, and they become much more computation- 
ally tractable. 



B. Tensor and Curvature Power Spectra 

With this newfound understanding, we explore the 
main spectral observables to see how they reflect the 
effective dimension and the main features of multi- 
field scenarios. We tackle these observables in rather 
quick succession, so we refer the interested reader to 
[El [m EU [23l [30] for more details about how these quan- 
tities are calculated. 

We start with the power spectra. The tensor power 
spectrum is unchanged by the presence of multiple fields 
and has the form [67] 



2^ 



(68) 



under the common convention for normalization of the 
spectrum. 

The density power spectrum is typically given in terms 
of the power spectrum of curvature perturbations, as the 
two are equivalent up to factors of 0(1) after inflation 
ends. The curvature perturbation IZ during inflation is 
related to the adiabatic density mode by [19] 



7^ 



(69) 



Using the mode sourcing equation for the adiabatic mode 
and quantizing the fields, one can find the curvature 



power spectrum. Using the notation in this paper, the 
curvature power spectrum at the end of inflation [23] can 
be rewritten as 



27r J 2e* 



-ns 



(70) 



where it is understood that the function T-jzs is evaluated 
at the end of inflation and where 



(iV) 



nN 

/ ' 



Dei 



dN 



TssiNi)dNi 



!2 p-/«.(lVIxx + Zxi)diVi 



(71) 



to lowest order in the slow-roll limit. We also made the 
replacement M + Z M±± + Z^^, which follows from 
the fact that the S(j)2 entropy mode can only be sourced 
by the other entropy modes. 



Examining Eq. (70), the curvature spectrum at 



the end of inflation equals the curvature spectrum at 
horizon exit plus corrections dictated by the function 
T7^5{^ The function Tss represents the evolution of 
the entropy- mode-related quantity taking into ac- 
count both the mode's intrinsic evolution and the effects 
of entropy mode sourcing. The other function, T-jzs^ is a 
measure of how much — sources the curvature mode — 

V 

or more liberally speaking, of how much 5(1)2 sources 5(l)i. 
The forms of the transfer functions follow directly from 
Eqs. (|48| and ([52|. We call Tns and Tss the multi-field 
transfer functions or simply the transfer functions since 
they are the multi-field generalization of the scalar trans- 
fer functions in two-field inflation [33l [661 • The difference 
is that there are multiple entropy modes in general multi- 
field inflation, so the transfer functions are necessarily 
vectors rather than scalar functions. 

Looking at Tss more closely, it involves the exponen- 
tial of the integral of a matrix, which can be estimated 
using the Magnus series expansion. According to the 
Magnus series expansion (see [68 and references therein) , 

if we take e"(^) = e" then the first 

three terms in the series expansion are 



We use the symbol Tt^^ rather than Nibbehnk and van Tent's 
symbol Up^ in order to mirror the transfer function notation 



used for two-field infiation p^[33ll66] . 
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/ (M + Z)^^^iiVi, 

N pNi 

[(M + Z)^^(7Vi), (M + Z)^^{N2)]dN2dNu 



no =- 



(72) 



^^=-tJ / / ([(M + Z)^i(iVi),[(M + Z)ii(iV2),(M + Z)ii(iV3)]] 

<5! Jat, Jat, Jat, 

Z)±±(iVi)]])dAr3diV2diVi, 



In, Jn, Jn, 

[(M + Z)^^{N3), [(M + Z)ii(Ar2), (M ■ 



where [A, B] = AB — B A is the matrix commutator of 
matrices A and B. Note that because M and Z are 
symmetric and anti-symmetric, respectively, the commu- 
tator [M^ Z] = 0, and hence the above expressions fur- 
ther simphfy, with the result being only commutators of 
each matrix with itself at different time points remain. 
The computational complexity of the Magnus series ex- 
pansion vividly illustrates how advantageous it is to use 
simplifications. For every mode degree of freedom that 
we can ignore, the calculation of the matrix exponen- 
tial greatly simplifies. Hence we strongly emphasize the 
importance of using our newfound result that M is a 
tridiagonal symmetric matrix in the slow-roll limit in the 
kinematical basis, of determining the effective dimension 
of a multi-field scenario, and of exploiting matrix sym- 
metry and anti-symmetry properties to simplify the com- 
mutators. 

Now the geometry and kinematics of the infiationary 
Lagrangian affect the curvature spectrum through the 
quantities H^^ 2e*, and Tt^^ (which in turn depends on 
The quantities and straightforwardly de- 
pend on the Lagrangian geometry since they are func- 
tions of K and |V*lnV|, respectively. The third quan- 
tity is more complicated. It depends on the turn rate 
of the background trajectory times the transfer function 
a vector function representing how much the S(j)2 
mode itself is sourced by the other d — 1 entropy modes. 
But fortunately, the behavior of Tt^^ can be inferred 
from our classical treatment of the mode sourcing equa- 
tions. For example, if the 62 basis vector is rapidly turn- 
ing into the 63 direction while ei turns at a constant rate, 
then 6(1)2 will be strongly sourced by ^^3, causing a boost 
in the amplitude of both transfer functions. As a second 
example, if the field trajectory rolls along and parallel to 
a ridge in the potential, then the S(j)2 niode will grow in 
amplitude, causing a boost in Tss and a greater than 
perhaps expected boost in Tt^^. As a third example, if 
a strong negative curvature Rs2 arises from the kinetic 
terms in the Lagrangian and dominates the dynamics of 
the 6(1)2 and ^^3 modes, both modes will decay, thereby 
reducing Tss and blunting the sourcing function Tt^^. 
This demonstrates that the effects of the geometry and 
kinematics of the Lagrangian on the curvature spectrum 
follow directly from all the connections we made in Sec- 



tion |III C| between the mode sourcing and the geometry 
and kinematics of infiation. Thus, we emphasize that one 
can understand how the Lagrangian translates into cos- 
mic observables by studying the mode sourcing in detail. 

Next, to determine how the effective number of fields 
is refiected in the curvature spectrum, we define a new 
unit vector 



e7^ 



and the scalar quantity 



-ns ' 



(73) 



(74) 



The unit vector necessarily lies in the {d — 1)- 
dimensional subspace spanned by the kinematical ba- 
sis vectors ej, 63, e^, where again * represents that a 
quantity is evaluated at horizon exit. If infiation has two 
effective fields, then = ej; however, if infiation has 
more than two effective fields, then 7^ Moreover, 
one plus the number of nonzero components of in the 
kinematical basis gives the effective number of fields. We 
use the above quantities to rewrite the power spectrum 
for general multi-field inflation as 



(75) 



Recall that in single-fleld inflation, T^s = 0, and the 
single-fleld curvature power spectrum is times larger 
than the tensor spectrum. That is, the tensor-to-scalar 
ratio tt, deflned by 



Vt'' 



produces the single-fleld consistency relation 
Vt = — 8nT, 



(76) 



(77) 



where the tensor spectral index ut = ~ — 2e*. In 

multi-fleld inflation, however, the ratio satisfles the upper 
bound UniEII: 



Tt = —Stitcos An < —Sut^ 



(78) 
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where 



tan A 



N 



(79) 



But beyond this, the curvature power spectrum cannot 
provide any further clues about the effective number of 
fields during multi-field inflation. The reason why is be- 
cause Eq. (75) is identical in form to the corresponding 



expression for two-field inflation [13 . To probe the num- 
ber of effective fields during inflation, we have to obtain 
information on the number of non-zero components of 
e-T^, which is not possible to do with only the tensor and 
curvature power spectra. 

As an aside, the mult i- field curvature power spectrum 
can also be given in terms of the 5N formalism. Un- 
der the 5N formalism, correlators of IZ can be written 
in terms of covariant derivatives of N ^ so the curvature 
power spectrum can be written as 



(80) 



where VA/" is the covariant derivative of the number of 



e- folds of inflation. By comparing Eqs. (70) and (80) and 
using that • e7^ = 0, it follows that 



1 



: (e^ + Tns^n) , 



(81) 



and therefore, the unit vector in the direction of A" is 

Gat = cos AjveJ + sin AatCt^. (82) 

These results generalize those found for two-field inflation 
in [69 and will be useful later when we calculate the non- 
Gaussianity arising from multi-field inflation. 

C. Isocurvature and Cross Spectra 

If there is more than one field present, then there 
will also be a relic spectrum of entropy modes and a 
cross spectrum between the density and entropy modes. 
Therefore, the detection of an entropy mode spectrum 
arising from inflation would indicate that inflation is 
multi-field. In lieu of the entropy modes, one usually 
works in terms of isocurvature modes <S, which can be 
defined in a gauge- invariant and dimensionless manner: 



~ p' p'' 



(83) 



Up to normalization factors, it can be shown that S de- 
pends on only the entropy mode 5(j)2- Here, we choose the 
normalization factor so that the isocurvature and curva- 
ture spectra have equal power at horizon crossing: 



S = 



(84) 



Under these above assumptions and using some prior re- 
sults from lini |23] , the isocurvature spectrum at the end 
of inflation is 



Vs = 



(85) 



where Tss is given by Eq. ( 71 ) and is calculated at the 
end of inflation. 

Now, unlike for the curvature modes, determining the 
isocurvature spectrum in the period after inflation ends is 
much more complicated. After inflation ends, the isocur- 
vature modes may decay further, affecting the amplitude 
of the isocurvature spectrum at the time of the CMB. But 
such post-inflationary processing is highly model depen- 
dent and depends on the dynamics of reheating, among 
other factors. So to keep our discussion as simple and 
as broadly applicable as possible, we focus only on the 
amplitude of the isocurvature modes at the end of infla- 
tion. Any post-inflationary modification of the isocurva- 
ture modes can always be tacked onto these results, and 
so these results can be construed as upper limits on the 
amplitude of isocurvature modes. 

Like in the previous section, how the geometry and 
kinematics of inflation affects the spectrum follows from 
our detailed discussion of the mode sourcing in Section 



III C| So let us proceed by considering how the number 
of fields is reflected in the isocurvature spectrum. Like 
for the other transfer function, we can break Tss into 
two parts: 



Tss 



IT55I 



(86) 



In the case of two-field inflation, = 63, whereas for 
inflation with three or more effective fields, 7^ ej. 
Using these two quantities, the multi-field isocurvature 
spectrum is 



2^ 



1 

2^ 



(87) 



Like for the curvature spectrum, the expression for the 
multi-field isocurvature spectrum has the same form as in 
the two-field case and therefore does not provide us any 
insight into the number of fields present during inflation, 
at least not to lowest-order in the slow-roll expansion. 

Also if inflation is multi-field, there will be a cross spec- 
trum between the curvature and isocurvature modes, rep- 
resenting the correlations between the two modes. Again 
combining results from [131 ES] , we write the cross spec- 
trum as 



Obviously, inflation is necessarily single-field if there 
is no cross spectrum, but let us consider if the cross- 
spectrum can provide any information into the number 
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of fields during multi-field inflation. Using Eqs. (73) and 
(86), the cross-spectrum can be written as 



Cns 



2^ 



1 

2^ 



(89) 



Now comparing the above result to [13], we see that the 
above result is identical to the two-field result with the 
exception of the term • ^s- This is the first instance of 
a spectral quantity whose expression differs from the two- 
field case and hence can potentially be used to gain more 
insight into the effective dimension of a multi-field sys- 
tem. If • = 1, as is automatically true for two-field 
inflation, then we recover the two-field result. However, 
if the two unit vectors are not parallel to each other, then 
there must be three or more fields present during infla- 
tion, and this should be reflected in the cross spectrum. 
We can therefore use the cross spectrum to devise a test 
that will distinguish inflation with three or more fields 
from inflation with two fields. In analogy to the tensor- 
to-scalar ratio, we define the cross-correlation ratio as 
[32] 



rc 



Cns 



which implies that 



rc = sin AA^e7^ • e^. 



Substituting Eq. (97) into Eq. (91), we obtain 



rc<x I 



SriT ' 



(90) 



(91) 



(92) 



where equality implies the presence of only two effective 
fields and inequality implies three or more fields. 

Based on these observations, it is therefore helpful to 
define two multi-field parameters 



Pi 

P2 



SriT ' 
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(93) 



The multi-field parameter f3i distinguishes single-field in- 
flation from multi-field inflation. It derives from the well- 



known single-field consistency relation in Eq. ( 77 ) . When 
/3i = 1, inflation is single- field; whereas, for multi- field 
inflation, < /3i < 1. The second multi- field parameter 
distinguishes among multi- field models, differentiating 
two-field models from models with three or more fields. 
When = 1, inflation is driven by two fields, whereas 
for models with three or more fields, < f32 < 1. Since 
we are not accounting for any possible post-inflationary 
decay of the isocurvature modes, in practice, P2 is more 
informative when it equals 1 (signaling two fields), rather 
than when it is less than 1 (signaling either more than 
two fields or post-inflationary decay of the isocurvature 
modes). These results are summarized in Fig. [sj 



Multi-Field Observables 



SriT 



Inflation 



0<Pi< 1 



Pi: 




Single-Field 



Multi-Field 



P2: 




0<p2< 1 



Three- or 
More-Field 



FIG. 3. Multi-field observables /3i and 132 indicate the effective 
number of fields during inflation. 



In sum, the tensor, curvature, isocurvature, and cross- 
spectra can in principle be used to distinguish among 
inflationary models driven by one, two, and three or more 
fields. However, the caveat is that all four spectra must 
be measured. 



D. Higher-Order Spectra 

Let us consider whether other spectral observables 
provide any clues about the number of fields present 
during inflation. Another group of important cosmic 
observables are the Fourier transforms of higher-order 
mode correlation functions. These higher-order spectra 
represent the non- Gaussian behavior of the perturba- 
tions. The two lowest-order correlation functions are 
known as the bispectrum and trispectrum, respectively. 
For standard multi-field inflation, the local form of 
these curvature spectra predominate}^ with the local 
bispectrum represented by the parameter /atl and the 
trispectrum by the parameters tatl and gNL- Using the 
6N formalism, the non-Gaussianity parameters can be 
expressed as [7lti73] 



^'^ However, when R 7^ 0, other forms of the bispectra and trispec- 
tra may also be important. We only give the local form here. 
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In [SS] , we calculated the above non-Gaussianity parame- 
ters for two-field infiation by operating V on the transfer 
function expression for V^AT in Eq. fell. The calcula- 



tion is similar for the case of general multi-field inflation. 
Repeating the steps we followed in [69 , we obtain a semi- 
analytic formula for the bispectrum parameter: 



-^/atl = ^cos^ ^N{nn - rir) + sin A at cos A at 



(e^Mei)* + sin Aat cos Aat a/— ' VT-j 
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The only difference between the above result and the re- 
sult for two-field infiation is that e2 has been replaced by 

To find the trispectrum parameters, we also follow a 
similar procedure to [69 . For tatl, we obtain 

n 2 
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sin"' A 
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fNL + \ cos^ AA^(n7^ - ut) 



Since 
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we can write tnl completely in terms of other spectral 
observables, which gives us a new consistency condition 
for general multi-field infiation: 
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which is only valid when infiation contains multiple fields. 
In the limit where /atl is detectably large (i.e., I/atlI ^ 
3), the above multi- field consistency condition reduces to 

2 
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In this limit, the value of tatl relative to /^^ is controlled 
solely by the ratio of to n^; the larger the sourcing of 
the curvature mode by the isocurvature modes, the more 

tnl approaches (I/atl) • In other words, only multi- 
field inflationary scenarios where the multi-field effects 
are very weak can produce tnl^ Inl- 

Lastly, for the trispectrum parameter Qnl^ we find 
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As written, the above result for multi- field inflation is 
identical in form to that in two-field infiation. 

Let us consider how the results for the three param- 
eters /atl, ^nLi and qnl reflect the effective number of 
fields. In single-field infiation, the expressions for the 
three parameters simplify nicely. In fact, the parameters 
fNL [70^, Tnl [74 , and qnl each give rise to a single-field 
consistency relation: 
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where ^ /atl represents the scale dependence of 

/atl and where we combined the single-field limits of Eq. 
(96) to obtain the last relation. Hence, a measurement 



of any one of the three non-Gaussian parameters can tell 
us whether infiation is single-field or multi-field. How- 
ever, these parameters provide no additional insight into 
the effective number of fields present during multi-field 
infiation. The reason why is that their expressions in 
multi-field infiation are identical to those in two-field in- 
flation after the replacement 63 67^, and hence they 
cannot differentiate models with two fields from those 
with three or more fields. 



V. CONCLUSIONS 

In this paper, we explored how the geometric and kine- 
matical features of multi-field inflationary Lagrangians 
with noncanonical kinetic terms are reflected in the cos- 
mic observables. We focused on how these geometric and 
kinematical features affect the interactions among modes 
and how this determines the effective number of active 
fields during infiation. 

We started by presenting the covariant background 
equation of motion for the fields and by delineating a 
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framework to parse the kinematics of the background 
fields. The kinematics of the background fields induce 
a basis called the kinematical basis and a matrix of turn 
rates, Z, which characterizes how quickly these basis vec- 
tors are rotating. This framework is important because 
working in the kinematical basis yields several advan- 
tages in elucidating the behavior of both the unperturbed 
and perturbed fields. We concluded our treatment of the 
background fields by discussing the slow-roll limit and 
certain subcases of limiting behavior. 

Next, we explored the equation of motion for the field 
perturbations in both the given and kinematical bases. 
In the combined super horizon and slow-roll limits, the 
equation of motion for the field perturbations depends 
only on the effective mass matrix M (a primarily geo- 
metrical tensor representing the covariant Hessian of the 
potential and the Riemann tensor of the field manifold) 
and the turn rate matrix Z (a kinematical tensor). We 
then studied the mode interactions one by one, filling 
an important gap in the literature. To do so, we first 
proved that most of the coefficients of the mass matrix 
M = V In V vanish in the kinematical basis. We used 
this result to derive a series of simplified mode sourc- 
ing equations for all the modes in the slow-roll limit. In 
particular, we found that for Lagrangians with canon- 
ical kinetic terms, the matrix sum M + Z is tridiago- 
nal in the kinematical basis, so the (5^^ entropy mode 
can only be sourced by the 5(j)n-i and S<pn-\-i entropy 
modes. For quadratic potentials with canonical kinetic 
terms, the equations simplify even further, in a way such 
that each mode S(j)n can be sourced only by ^^n+i but 
only when the basis vector is turning into the direction 
of en+i- We therefore argued that the mode interactions 
in a general inffation model can be divided into features 
shared in common with canonical quadratic models and 
features that arise from differences from this fundamental 
Lagrangian, and we advocated this approach as way to 
gain greater insight into how differences in Lagrangians 
translate into differences in the cosmic observables. 

In parallel, we discussed how the geometric and kine- 
matical features of multi-field infiationary Lagrangians 
affect the mode sourcing. We started by considering the 
case where the S(j)n mode is unsourced, and we discussed 
how the concavity of the potential and the curvature of 
the field manifold determine that mode's instrinic evolu- 
tion rate. We showed that there are up to d mode-related 
quantities in infiation with d scalar fields that may be 
conserved, thereby generalizing the single-field conserva- 
tion law for the adiabatic mode. Next, we discussed how 
the sourcing effects for the S(j)ri mode depend two types of 
geometrical terms and one type of kinematical term. The 
geometrical terms involve off-diagonal terms in both the 
covariant Hessian of the potential and in the Riemann 
tensor of the field metric contracted with two instances 
of ei and modulated by e, and we interpreted these terms 
geometrically. The kinematical terms are simply the turn 
rates of into the e^_i and e^+i directions. Although 
the turn rates can be expressed in terms of geometri- 



cal features of the Lagrangian, these expressions are in 
general complicated and opaque, so it is usually easier 
to treat these terms as kinematical quantities. We also 
gave several examples of how inferences about the mode 
sourcing can be made by determining the geometric and 
kinematical features of a Lagrangian. 

Determining how these geometric and kinematical fea- 
tures are refiected in the cosmic observables mostly boils 
down to understanding how the modes interact in the 
classical picture. Therefore, it is essential to study the 
interactions among modes as we demonstrated. With 
this in mind, we then focused on how the Lagrangian ge- 
ometry and kinematics determines the effective number 
of fields and how this number is refiected in the power 
spectra, bispectrum, and trispectrum. We pointed out 
that the effective number of fields needed to describe the 
background and perturbed solutions do not necessarily 
coincide, and we gave a simple method to determine the 
effective dimension of a multi-field system in the slow-roll 
limit. Next, we presented known formulas for the power 
spectra, discussed how the simplifications presented in 
this paper can greatly ease the calculation of the spec- 
tra, and then generalized the two-field expressions for the 
local non-Gaussianity parameters to multi-field infiation. 
We extended our two-field consistency relation [69] relat- 
ing several of the cosmic observables to multi-field infia- 
tion, and we discovered a multi-field observable involving 
the cross spectrum that can potentially distinguish two- 
field models from models with three or more effective 
fields. However, being able to use the cross spectrum for 
this purpose requires knowledge of the post-inflationary 
history of the entropy modes. In the future, a better 
understanding of this post-inflationary evolution is es- 
sential in order to gain better insights into the spectral 
observables involving the isocurvature modes. 

Stepping back and looking at the big picture, since 
more sensitive measurements of the spectral observables 
plus new spectral observables will reveal further clues 
into the nature of inflation, we must be posed to extract 
phenomenological information from these measurements. 
However, it is impractical to test the myriad inflationary 
models one by one against these measurements. Instead, 
we advocate studying classes of geometric and kinemat- 
ical behaviors that arise from inflationary Lagrangians 
and determining how these behaviors affect the cosmic 
observables. This understanding can then be reversed 
and used to determine how constraints on spectral ob- 
servables in turn constrain inflationary models, so that 
we can eventually put constraints on the inflationary La- 
grangian that described our early Universe. The work 
presented in this paper represents a key step forward to- 
wards this ultimate goal. 
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